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A simple model for the peeling process of pressure-sensitive adhesive tape is presented. 
The model consists of linear springs and dashpots and can be solved analytically. Based on 
»-G , the modeling, the curved profile of the peeling tape is spontaneously determined in terms of 

viscoelastic properties of adhesives. Using this model, two experimental results are discussed: 
C^ , critical peel rates in the peel force and the peel rate dependence of the detachment process 

of adhesive from the substrate. 
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►-j^ ' 1. Introduction 

^2 ■ The behaviors of the peeling process of pressure-sensitive adhesive tape are very complex 

t^^ ■ and interesting: (i) The peel force curve against the peel rate sometimes has two or three 

^^ ' local maximums (Fig. 1).^'^ This fact implies that we will observe self-excited oscillation of 

f^ I the peel force, if the adhesive tape is peeled at a rate in the negative slope in Fig. 1.^'^ (ii) 



C^ 



We observe fingerlike instabilities at the peel line, which is the interface between the adhesive 

and the rigid substrate (Fig. 2), and the front of failure waves along the peel line.^ (iii) In 

'^ ' some adhesives, we observe caves inside the adhesive.^ The appearance of the caves indicates 

C ■ 
Q I the flow of the adhesive; the adhesives are polymer melts with viscoelasticity. 

In this paper, in order to describe the peeling system with viscoelastic adhesives, we 

present a simple model made of linear strings and dashpots and solve it analytically. A few 

^ • attempts have already been made to incorporate the viscoelastic properties.^' ^ The advantage 

C^ ' 

of our modeling is on the spontaneous determination of the curved profile of the peeling tape 

with the viscoelastic adhesives. 

Taking the advantage, we shall discuss the followings: (i) "critical peel rate" at which the 

peel force curve takes a local maximum and (ii) the positive dependence of the peel force 

on the peel rate, which is not trivial. Actually, it is not straightforward to derive positive 

dependence, and we approach this difficulty by considering more microscopic picture at the 
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interface between the adhesive and the substrate in § 4. 

2. Model 

2.1 Frame of the model 

We shall consider the system in two dimensions of the x-y plane in Fig. 2. In practice, 
there is z-direction dependence of the system with fingerlike instabilities along the peel line. 
However, with small amplitude of the fingering, this assumption will be justified. 

The adhesive is represented by a set of linear springs (A) with spring constant k/N where 
A^ is the number of the springs per unit length along the x axis. These springs are numbered 
along the x axis as indicated in Fig. 3; k corresponds to the modulus of elasticity of the 
adhesive Ea and will be represented as /c = EaOJJy^ for the adhesive with the width lo and 
initial thickness y^. These springs can be elongated only along the y axis. We shall call this 
part of the model just "element"; the element will be replaced by another element, such as 
a series connection of linear spring and dashpot (the so-called Maxwell element) or various 
linear combinations of linear springs and dashpots. 

Secondly, corresponding to the backing tape shown in Fig. 2, we consider flexible joints 
which connect adjoining springs (A) and are represented by linear springs (B) with spring 
constant NK^ as shown in Fig. 3. This modeling assumes that the backing tape is a thin 
elastic object without bending; K corresponds to the modulus of elasticity of the backing 
tape i?b and will be represented as i^ = Ei,ijjb for the backing tape with the width uj and 
thickness b. We admit that it is more realistic to describe the backing tape as a beam with 
flexural rigidity. ^^ However, if we adopt the bending of the beam, the equation system becomes 
much more complex, as shown in Appendix A, and the results are not essentially different 
from the results with the flexible joint of the linear springs (B). Therefore, we adopt this 
simpler model of the linear springs (B) as the flexible joint in the main part of the following 
discussion. In terms of the rigid substrate, it is represented by a horizontal line (C) (the x 
axis) as in Fig. 2. It is assumed that all elements and springs (B) are massless. 

We impose on the elements a condition that the bottom of each element is detached from 
the line (C) when the force applied to the element becomes some critical value (Tc/N . Imposing 
this condition on the element corresponds to the situation where the adhesive is detached from 
the substrate when the normal stress of the adhesive at the interface reaches the critical value 
Cc/w. Until § 4, we assume that the quantity Uc is constant; the possibility of the peel rate 
dependence of ac is discussed in § 4. As a natural boundary condition, we assume there is no 
force at x = oo. 

Under these assumptions, let us now consider the following situation: the linear springs 
(A) and (B) are settled in the range [—00,00] and there is no force acting on the springs in 
the initial stage. At t = —00, we begin to lift up quasi-statically one spring (B) at the left end 
at X = — 00 ; since we move it quasi-statically, the force balance holds for every element at 
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every moment. We suppose that the force apphed to the element at x = reaches the critical 
value Gc/N at t = (we call such an element the "edge" of the elements hereafter because it 
is the boundary between the elements having left the substrate and the elements sticking to 
the substrate). 

Under such situation, let us now try to obtain the shape of the curve drawn by the top 
ends of the elements (we call this curve "flexible joint's curve" for brevity). First, we consider 
the force balance equation at the i-th element (i > 0): 

-^fi - NK{h - h+i) + iVK(/i_i - fi) = 0, (1) 

where fi is the y coordinate of the top end of the i-th. element. From the conditions at i = 
and i = oo, we obtain the following relations: 

^/o = I (2) 

/oo = 0. (3) 

Taking the limit A^ ^ oo in these relations and rewriting the coordinate fi by f{x) through 
the relation x = -^, we obtain the equation for the flexible joint's curve and its boundary 
conditions from (l)-(3): 

Kf"{x) = kf{x) (4) 

Kf'iO) = a, (5) 

/(oo) = (6) 

where the prime over the function / means its derivative with respect to x. Solving eq. (4) 
under the boundary conditions (5) and (6), we have 

K 

where a = y-^; a characterizes the range where the elements are stretched (we call this 
length the "stretching length" for brevity). 

From the force balance relation at the 0th element, the peel force F applied to the spring 
(B) at X = is represented by the flexible joint's curve /, as follows, 

F = ^f^ + NKUo-fi). (8) 

Then, after taking the limit A^ — > oo, we have 

F = -Kf'{0) (9) 



fix) = ^e-^/-, (7) 



which implies for this model that 



F = aoc. (10) 
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It is noted that the flexible joint's curve obtained in our modehng with more complex 
elements comprised of linear combination of springs and dashpots is also of the exponential 
form as 

/(x) = /(0)e--/^ (11) 

where /(O) and (3 are some positive constants; the proof is given in Appendix B. From the 
boundary condition (5) which is satisfied with more complex elements, /(O) is turned out to 
be /(O) = [5'^cjf./ K. Then, with the general relationship (9) which is also satisfied with more 
complex elements, the peel force always has the form of 

F = (3ac. (12) 

This relation means that the peel force is given by the multiplication of the critical stress Cc 
and the stretching length (3. In the present modeling, this stretching length (3 is spontaneously 
determined in terms of the elastic moduluses of the adhesive, k, and of the flexible joint, K. 
Therefore, when we keep both the elastic modulus of the flexible joint K and the critical stress 
Gc constant, the peel force monotonically decreases as the adhesive gets stiffer with larger k. 
This relation becomes important in understanding the behavior with viscoelastic elements 
showing V dependence of /c, as discussed in the following. 
2.2 Case of Maxwell elements 

As mentioned in § 1, since the adhesive is a polymer melt with visco-elasticity, it is more 
plausible to describe the adhesive in terms of Maxwell element of serial connection of spring 
and dashpot with spring constant, /c, and viscosity, r/. For the serial connection of spring and 
dashpot, the force acting on the element is given as kyi = r}y2, where yi and y2 represent the 
displacements of the spring and dashpot, and the dot over y2 rneans its time derivative; the 
total displacement y is given as y = yi + y2- The general solution of this system gives the 
applied force as k j_ e ~y{t')dt' , where r = r]/k is the relaxation time of the Maxwell 
element. 

Suppose that the spiring (B) at the left end is lifted up with constant speed v a.t t = —oo 
and the system is in the steady state, then the y coordinate of the i-th element sticking to 
the substrate will be expressed as f{i/N — vt). If the edge of the elements reaches x = at 
t = 0, we get the following equation of f{x) for the limit of A^ ^ oo from the force balance at 
the i-th element (i > 0): 

Kf"{x) = k e^'/^f{x - vt')dt', (13) 

J —oo 

where x corresponds to i/N. It is useful to note that the visco-elastic properties of the element 
appears only in the right hand side of the above equation. The boundary conditions are the 
same as eqs. (5) and (6). 
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Suppose that the solution of this equation has an exponential form of eq. (11) and substi- 
tute it into eq. (13), then we get an algebraic equation for (3: 



2 

p^-—p-a^ = 0. (14) 

VT 



where the condition /? > is set by the boundary condition (6), and hence 



/?= — + -\/(— )2 + 4a2. (15) 

Using the boundary condition (5), we get f{x) = {I3'^ac/K)e~^'l^ and we obtain the peel force 
represented as eq. (12) from the relationship (9). We have plotted in Fig. 5 the peel force F 
against the peel rate v for r = 1, /c = 1, i^ = 1, and Uc = 1. 

For the case of Maxwell element with viscoelasticity, the exponential deformation brings 
an additional feature that the whole elements behave in a similar way; i.e. all elements behave 
as a viscous liquid at low peel rate and as an elastic solid at high rate, as shown in the following 
expression of the applied force to each element, 

k f e*'/V(x - vt')dt' = -^^kf{x) = ,^^nfix) (16) 

j_oo 1 + vra 1 + vra 

In this sense, this system of peeling adhesive tape is essentially different from the system of 

a fracture of adhesive studied by de Gennes,^^ where the behavior of an viscoelastic adhesive 

depends on the position in the system. As the consequence of this behavior at the limiting of 

V ^>- 0, the stretching length (3 expressed by eq. (15) and consequently the peel force F diverges 

at t> = 0, as shown in Fig. 5. The divergence of the peel force at f = is obviously unrealistic. 

In the next section, we introduce more realistic model with a combination of springs and 

dashpots to describe the adhesive of polymer melt with slight cross-links. 

We finally comment on the stability of this system. In the present analysis, the expression 

of /(x — vt) assumes the instantaneous adjustment to the stationary state and hence assumes 

the stability. From the stability analysis of an analogous situation,^ we can put one comment 

to justify the assumption. If the whole system is hard enough to transfer the change in applied 

load to the elements in a shorter time interval in comparison with the relaxation time of the 

elements, the faster rate of detachment of the elements immediately brings smaller applied 

force on the elements under the condition of constant speed of lifting so that the detachment 

rate becomes slower, and vice versa. That is, for the hard system, the peeling edge tends to 

move with constant speed equal to the lifting velocity. 

3. Explanation of Two Critical Peel Rates in Terms of the Model 

As mentioned above, the adhesive is a polymer melt with small portions of cross-linking. 
The visco-elastic properties of the response to the deformation can be represented by two 
characteristic times; they are the times characterizing the transitions among the glass-like 
response for fast deformations, the rubber-like responses with entanglements behaving as 
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temporal cross-links and with true cross-linking for intermediate and slow rates of deformation, 
respectively, which are denoted by Tg and r^ (r^ < r^), respectively.^ 

Thus we should construct the element of the model so as to have two relaxation times, as 
shown in Fig. 6, which consists of one linear spring with spring constant k and two different 
Maxwell elements which are specified by the elastic constants, kr and kg, and viscous constants, 
T]r and r]g, respectively; these quantities are related to the relaxation times Tr and Tg as: 
Tr = r]r/kr and Tg = rjg/kg. We call each part of the element "spring (I)", "Maxwell element 
(II)" and "Maxell element (III)," respectively. Note that spring (I) represents the effect of the 
cross-linked polymer within the adhesive. 

The equation for f{x) is 

/■O /-O 

/"(x) = a-^f{x) + ar-'^ e''/^^f{x-vt')dt' + ag-^ e^' '^^ f{x - vt')dt' (17) 



where a = \J Kjk, Or = y^K/kj., and Og = y/K/kg. The boundary conditions are the same 
as eqs. (5) and (6). 

The peel force is given as eq. (12). The stretching length /3 is the positive root of the 
following algebraic equation 

/?"^ = a"^ + a-'^VTr/if3 + VTr) + a~'^VTg/{l3 + VTg). (18) 

Plotting the peel force F against the logarithm of the peel rate v, as shown in Fig. 7, we can 
find two critical peel rates labelled by vi and V2- These two peel rates vi and V2 corresponds 
to the rates at which the elastic modulus of the element gradually changes from k to k + kj., 
and from k + kr to k + kr + kg, respectively, with increasing peel rate. 

At a very low rate. Maxwell elements (II) and (III) do not contribute to the elastic modulus 
of the element because the elements have already been relaxed. Thus the elastic modulus 
of the element turns out to be k only, which is contributed by the spring (I). From the 
result around (7), the "stretching length" in that case is given by ^J Kjk. Then, the time 
required for elongating each element until its detachment from the substrate is estimated to 



be {l/v)\/ K/k, where v is the peel rate (we call the time "stretching time"). As the peel 
rate increases, the stretching time {l/v)yK/k gets closer to the relaxation time of Maxwell 
element (II), Tr, and hence the elastic modulus of the element gradually changes from k to 
k + kr- Thus the peel rate vi is estimated as, 

"■ = tM '■''^ 

By the same logic, we can also estimate the peel rate, V2 as. 



In Fig. 7; we can actually see that the evaluated f i and V2 correspond to the boundaries of 
the ranges. 
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As shown in Fig. 7, in our model with constant critical stress for the detachment cJc, 
the peel force curve is a decreasing function of the peel rate and does not reproduce the 
experimental results of the positive slope. If ac is an increasing function of the peel rate, the 
peel force curve will then recover the positive dependence on peel rate with small regions 
having negative slopes around the two critical rates vi and V2. Then, the peel force curve will 
become very similar to the experimental results schematically shown in Fig. 1. 

Of course, we cannot readily say the reason for the dependence of the critical normal 
stress (Tc on the peel rate (we are giving one possibility of the peel rate dependence of the 
critical stress in the next section), but it is probably obvious that the appearance of the 
two critical rates in the experiment is due to the changes in elastic modulus of the adhesive. 
Actually, by a lot of investigators, this relationship between the change in elastic modulus of 
the adhesive and the critical peel rates has been postulated.^ We have succeeded in showing 
the relation between the critical peel rates and the relaxation times of the adhesive in terms 
of the quantities characterizing the viscoelastic properties of the adhesives. 

4. On the Peel Rate Dependence of the Critical Normal Stress <Tc 

As discussed in the previous section, if we calculate the peel force in our model with 
constant critical stress ac for detachment, we inevitably observe negative slope of the peel 
force, which is completely different from the experimental results even in the qualitative 
sense. Thus we must discuss the reason for the discrepancy in the slope in more detail. 

The element of our model is a reduced one representing some part of the adhesive with 
finite size, so that the element can be considered to consist of more microscopic parts, such as 
polymer chains. So we can consider the bottom ends of the elements as the assembly of more 
microscopic parts as illustrated in Fig. 8 (we call them "small springs"). 

Since these small springs are microscopic objects such as polymer chains, it is reasonable to 
assume that their processes of detachment and attachment are stochastic; each small spring is 
stochastically detached and attached to the substrate with some probabilities which generally 
depend on the force applied to the small spring. 

We expect that the small springs are detached from the substrate more often as the applied 

force gets larger. We may also expect that the attachment process of each small spring depends 

on the applied force. But the dependence will be weaker, because the detached small springs 

will stay close to the substrate and will attach to the substrate almost equally for finite applied 

force. Thus we shall write down the time evolution equation for the amount A of small springs 

being attached to the substrate (which is a normalized quantity) as: 

dA{t) ^ A{t) l-A{t) 

dt T^{a{t),A{t)) r+ ' ^^ 

^ G{a{t),A{t)) (22) 
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T.{a,A) = roexp[— ^-^] (23) 

r+ = TO (24) 

where a is the tensile force acting on the whole element; a/ A is the force acting on each small 
spring, ^ is a constant having the dimensions of length, which corresponds to the width of the 
potential hole representing the attachment of small springs to the substrate, tq is a constant 
with the dimension of time, [/ is a quantity representing the depth of the potential minimum, 
kB is the Boltzmann constant, and T is the temperature of the system. When the quantity A 
vanishes the element leaves from the substrate. 

We first examine the stationary case, i.e., where a is independent of time. Solving the 
equation G{a,A) = for A and remembering that the sign of dG{a,A)/dA determines the 
linear stability of the solution, we readily find that the eqution system has its reasonable 
solution only when a is smaller than some critical value ag determined by the two equations 
G{as,A) = and dG{as, A)/dA = 0; namely, ag is given by the following equation, 

as + Hl^) = U-kBT. (25) 

If a is smaller than ag, A finally reaches a finite value in the steady state and the element is 
kept sticking to the substrate. If a is larger than as, on the other hand, A finally vanishes and 
the element is detached from the substrate. 

The behavior with time-dependent a is different from the stationary case. In our model, 
since the flexible joint's curve always has an exponential form, each element is subjected to a 
varying exponentially with time; i.e. a{t) = e^ *, in which v' is a parameter of a positive value. 
In such situation, a at the detaching moment, ac, depends on the rate v' because of the finite 
relaxation time in the detachment process; A does not vanish immediately. We plot cJc against 
the logarithm of v' in Fig. 9, where we have taken as the initial condition A{t = 0) = Aq 
where Aq is the stationary value of A for cr = because it is natural to assume that there 
is no force action on each element before the peeling process starting (the explict value is 
^0 = 0.999547). Owing to the positive dependence on v' of cJc, the positive slope will appear 
in the peel force curve as we usually observe experimentally. 

5. Discussion and Conclusion 

We presented a simple model to describe the peeling system, which is only made of linear 
springs and dashpots and solved the system analytically. With this modeling, the curved 
shape of the tape during the peeling is spontaneously determined in terms of the viscoelastic 
properties of the adhesive tape. The advantage of this characteristic is utilized in the discussion 
of the two critical peel rates being related with the viscoelastic properties of the adhesives. 

We discussed in § 4 the possibility of the peel rate dependence of the critical stress for 
detachment; if ac is constant we can only predict negative slope of the peel force which is in 
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contradiction with experimental results. In order to explain the dependence, we considered 
the microscopic picture of the element consisting of small springs repeating attachment and 
detachment to the substrate. We admit that there may be other possibilities for the peel rate 
dependence of the "effective" critical stress, e.g, the z dependence of the peel line with large 
amplitude of the fingering, etc. We have to carefully consider those possibilities theoretically 
and experimentally. 
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Appendix A: Case of Bending Beam as the Backing Tape 

We consider the backing tape as a beam with flexural rigidity Ki, = E^Iyy, where E^ 
represents the modulus of elasticity of the backing tape and lyy is the moment of the inertia 
of area and is expressed as lyy = w5^/ 12 for the backing tape with the width uj and thickness 
5. The differential equation and its solution for the case of elastic adhesive has been obtained 
by Bikerman^'^ in the following way. Firstly, the differential equation of the courved profile of 
the beam f{x) is given as,"^^ 

^'^d^ = -kf{x), (Al) 

where k has the same meaning as in our modeling of the adhesive with springs (A): k = EaU>/yo 
for the adhesive with the modulus of elasticity Ea, the width lo and initial thickness yo. The 
boundary conditions are 

Kbf""{0) = -a,, (A.2) 

/(oo) = (A-3) 

and additionally from the moment of force and its balance at x = 

/"(O) = 0, (A.4) 

F = Ktf"'{0). (A5) 

Solving the differential equation (A-1) under the boundary conditions (A-2)-(A-4), we have 

/(x) = ^e-"/"''cos(xK), (A6) 

where at = {AKh/k)^''^. The peel force F is then given from eq. (A-5) as, 

F = f cT, (A.7) 
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In the case of Maxwell elements having coefficients k and ry, we can construct the equation 
system and solve it in a similar way. The differential equation for f{x) is 



dx^ ^ " J-oo 



l{x) = - I e*/^/(x - vt)dt (A-9) 

The boundary conditions and the relationship between F and / are the same as (A-2)-(A-5). 
Supposing that / has the form f{x) = Ae~^'^ and substituting this form in (A-8), we 
then have an algebraic equation for /?, 

4(5^ + ^f3 + at = (A-10) 

VT 

Note that the real part of the root of this equation must be positive to satisfy the condition 
(A-3). We can find two such solutions; we denote them as P^ = uj + ij and P2 = uj — ij, 
where both uj and 7 are positive. Then, we can write / as 

/(x) = ^e""^^ cos(72: + 6), (A-11) 

where 9 = cos'^ J^ [0 < < tt) and 

A = -4 ^— (A-12) 

2cJ7(w2 + 72) Kb ^ ' 

which have been determined from (A-2) and (A-4). We can calculate the peel force from (A-5), 

(A-11) and (A-12), 

F = i-CTe (A-13) 

We plotted the peel force against the peel rate v in Fig. 10 for K^ = 1.25, k = l,r = 1 
and (Tc = 1. It should be noted that this result is essentially the same as the result of our 
modeling with the flexible joints of springs (B) shown in Fig. 5 in the following sense; (i) the 
peel force monotonically decreases with the peel rate; (ii) the peel force is saturated around 
some characteristic peel rate, which is estimated as (l/r)(4i^j,/A:)^'^. (In the case of springs 
(B) it is estimated as {l/T){K/ky/'^ .) 

Appendix B: Proof of the Form f{x) = /(0)e~^/^ for the Flexible Joints of 
Linear Spring (B) 

In our model with the linear springs (B) as the flexible joint , whatever the linear 
combination of linear springs and dashpots, the equation for f{x) necessarily takes the form 
as: 

M 

f"{x)=af{x) + bf{x-vt)_+^Ci //^^f{x-vt')dt'. (B-1) 

Here, a, 6, and Cj(z = 1, ..., M) are some non- negative constants and at least one of them must 
be nonzero, Tj is the i-th relaxation time of the system, and the dot over / stands for the 
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differentiation of / with respect to time. For a parallel combination of spring, dashpot and 
Maxell elements, there is a direct correspondence between the elements and the terms in the 
r.h.s. of (B-1); the first term corresponds to the spring, the second term dashpot, and the 
third term Maxell elements. For more general combination, the expression of (B-1) still holds 
though there is no such correspondence. 

Assume that the solution of (B-1) takes the exponential form of eq. (11) and substitute 
this form in (B-1), then we get an algebraic equation for /?: 

M 

VTi 



p-^ = a + bvr^ + y^Ci-—^. (B-2) 

^ P + VTi 



We prove that the algebraic eq. (B-2) has one and only one positive solution. Let us first 
introduce a variable ^ = 1/(3 and consider a function of ^ defined on the interval [0, oo], 

M 

^-^ 1 + VTiE 

We can readily find that g"{C) > for all non-negative ^^ i.e., 5 is a convex function of ^ on 
that interval. On the other hand, for a = 0, we see that (7(0) = and g'{^) < because at 
least one of b and Ci must be positive. For a > 0, we see that g{0) < 0. Noting g{oo) > and 
remembering that (7 is a convex function with respect to ^, we readily see that the equation 
5(0 = has one and only one positive solution and accordingly so does the eq. (B-2) for /?; 
we thus see the existence of a solution taking the form (11). 
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Logv 



Fig. 1. Schematic picture of the peel force against peel rate observed experimentally. The solid line 
indicates the peel force when the system is in the steady state, while the dashed line indicates the 
average value of the peel force; oscillation never stops in that rate region. 




Fig. 2. Schematic illustration of the experimental system. 
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Fig. 3. Schematic picuture of the prototype of the model with spring elements. 
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Fig. 4. Schematic picuture of the model in the case of Maxwell element. 



13/17 



J. Phys. Soc. Jpn. 



Full Paper 



Fig. 5. The peel force against the peel rate with constant ad—l) obtained in the case of 
linear springs (B) with A' = 1; where the element of the model is a Maxwell element 
with T = 1 and k = 1 
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Fig. 6. Schematic picuture of the element representing adhesive. The element has two 
relaxation times and the elasticity due to cross-links within the adhesive. 
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Fig. 7. The plot of peel force against the logarithm of the peel rate calculated from eqs. (12) 
and (18). Parameters are Tg = 0.00001 , t^ = 1 , iC = 1 , fc = 1 , fc^ = 1, (Tc = 1, and 
kg = 1. The gray line indicates vi and V2 evaluated from eqs. (19) and (20). 
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Fig. 8. Enlarged schematic picture of the bottom of the element. We shall consider the 
bottom of the element as an object consisting of small springs which are supposed to 
be detached and attached to the substrate stochastically. 
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Fig. 9. Plot of (Tc against logw' calculated from eq. (21) with a{t) 
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The parameter 



values in the calculation are r = 10, fcsT = 1, £ = 1, and U = 7.7. For these parameter 
values. CTc = 5.0734. 
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Fig. 10. The peel force against the peel rate with constant (Tc(=1) obtained in the case of 
bending beam with Kb = 1.25; where the element of the model is a Maxwell element 
with r = 1 and fc = 1, corresponding to Fig. 5 for the flexible joints of linear springs 
(B). 
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